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The CW microcavity optical parametric oscillator (OPO) state is investigated using a theoretical 
treatment which includes the contributions of the pump, signal and idler populations to the renor- 
malisation of the polariton energies (the 'blue-shift'). The theory predicts the pumping conditions 
under which the OPO switches on, showing that a pump angle >~ 10° is required, but there is no 
particular significance to the 'magic-angle' where pump, signal and idler are all on resonance. The 
signal and idler renormalisation contributions also causes the signal emission to be driven towards 
the normal direction as the pump power increases above threshold. 

The behaviour of a resonantly pumped semiconductor microcavity can be best described using the terminology of 
the optical parametric oscillator (OPO): with pumping above threshold at finite angle, a coherent signal appears, 
close to the normal direction, accompanied by an idler on the other side of the pumpii^ However, this is an unusual 
OPO, as the nonlinearity is x^i rather than the typical x^- One consequence is that there is a self-interaction, the 
'blue-shift', which can be thought of as a renormalisation of the polariton energies. This renormalisation obviously 
depends on the pump, signal and idler populations. However, although the term has been discussed in previous 
theoretical work, its effect on the OPO state have been ignored, or only the pump contribution has been considered. 
The purpose of the present paper is to show that, when the renormalisation is properly treated, a number of puzzling 
aspects of the experimental results can be explained. 

Recent experiments^ have shown that the CW-OPO can be switched on for a range of pump angles and energies, 
from 10'^ to at least 24°, with relatively small changes in threshold. This CW situation contrasts with the ultra- fast 
behaviour, where the the response is strongly peaked about a 'magic- angle', of ~ 16°, when the pump, signal and 
idler all lie on the polariton dispersion^i^ By investigating the stability of a state in which only the pump polariton is 
occupied, the theoretical treatment with the self-interaction provides a prediction of values of the pump angle, power 
and energy at which the OPO can occur. In fact, the OPO is not the only instability of the system: a simple bistable 
behaviour can also be obtained, as observed experimentally by Baas et aim for normal incidence pumping. The theory 
predicts that for small pump angles, only bistability occurs, while for angles greater than ^ 10° the OPO switches on. 
Furthermore, the OPO threshold is found to vary smoothly with pump angle and energy, giving no special significance 
to the magic angle. 

The other aspect of the experimental behaviour which is addressed here concerns the direction of the emerging 
signal, which is always found to be within a few degrees of the surface normal, whatever the angle of the pump. This 
is surprising because, if the signal angle were simply determined by the requirement that the mismatch of the signal 
and idler energies from the polariton dispersion is minimised, a wider range of angles would be obtained, depending 
on the pump energy and angle. The present theory predicts that this will in fact happen, but only very close to the 
threshold: at higher pump powers the renormalisation contributions from the signal and idler switch off the high angle 
OPO and pull the small angle states into resonance. In a theory without these contributions, the pump population is 
pinned to its threshold value, so nothing can change at higher pump powers. 

The first theoretical models of the microcavity OPO^ showed that the polariton energy renormalisation arises 
naturally from a description in terms of exciton-exciton scattering^ or a nonlinearityii Although the renormalisation 
was derived, its effects were not included in the solution for the OPO states. The renormalisation of the pump polariton 
state was shown by Baas et a& to explain the observed bistability, but this treatment was not extended to consider the 
OPO state. Gippius et aim discuss the bistability and OPO instability, making comparisons with a numerical model, 
including a continuum of modes, from which emerges the property of emission close to the normal. The present work is 
an analytic treatment which provides a physical explanation for the numerical and experimental results, by including 
the contributions to the renormalisation from the signal and idler populations in the description of the OPO. 

The treatment assumes that the pump beam is spatially uniform, and the pump, signal and idler fields arc simple 
plane waves. Although this is not fully realistic, as all experiments use a finite spot, it should give reasonable results 
for illumination with gradual spatial intensity variations. As will be discussed below, some differences between the 
theoretical predictions and experiment are probably due to the effects of spatially inhomogeneous pumping. It should 
also be noted that even for uniform pumping, it is possible that there are solutions with spatial structure to the 
internal fields. 

The theoretical expressions derived in this paper are fully two-dimensional, valid for any in-plane wave-vectors 
consistent with momentum conservation. However, for simplicity, the results presented in the figures assume the 
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FIG. 1: The pump-only state for 9p — 16°, uip = —2-y. (a) Detuning of the pump relative to the polariton energy, including 
the renormalisation contribution, (b) S-curve showing the relationship between the pump polariton population \p\^ and the 
external pump power Ip. (c) Values of \p\^ (shaded region) for which the pump-only state is unstable, as a function of the 
signal angle, 6s. Ii and I2 are the pump powers corresponding to the lowest and highest unstable values of |p|^ for any 6s: the 
OPO switches on when I2 > Ii- 



pump, signal and idler lie on a line passing through A; = in wave-vector space, and the excitations considered in 
the stability treatments are similarly confined. Calculations have also been made without this restriction, and stable 
OPO states which do not quite fall on such a radial line can be found. Similarly, extending the stability treatment to 
consider all excitations tends to reduce slightly the stabilty ranges of the solutions described. 

The remainder of the paper is developed in the following way: In Section HI the state with only pump polariton 
mode occupied is described, and its stability is investigated, leading to a criteron for determining the conditions under 
which the OPO switches on. In Section ^1 the model of the OPO state with finite amplitude is derived, showing 
how the emission is pulled towards the normal direction. Finally, in Sections IIIII and llVI the OPO state stability is 
investigated, and the various predictions of the model are discusssed. 



I. THE PUMP STATE 



This section describes a state in which only a single, pumped polariton mode is populated. The conditions under 
which this state becomes unstable are derived, and criteria are obtained for determining whether a simple bistability 
results, or the OPO switches on. 

The pump mode is assumed to be a plane wave cfip = (p/jXpp) expi(kpX — t^pt), where the amplitude, p, satisfiesii 

1 |C I 

\Ap\ 

Here, uJo{kp) is the lower branch dispersion at the pump wave- vector, 7^ its width, \ Cp \ and \Xp \ the cavity and exciton 
amplitudes (Hopficld factors), and fp the external pump amplitude. The non-linear term, represents the self 

interaction of the pump polaritons. By appropriate scaling of the fields^ the coefficient k can be made equal to unity, 
which forthwith is taken to be the case. The use of this single-branch form assumes that the polariton coupling is 
sufficiently large that the two branches are not mixed by the non-linear terms. 

The solution of the cubic Eq^has been discussed in ReflQ): when the detuning ujp — ujo{kp) > VS"fp, the cubic has 
two turning points, resulting in a characteristic S-curve, as shown in Fig^b)ii. For pump-powers falling between 
the two turning points, the system is bistable (the branch with negative slope is unstable). Fig^a) shows how 
this bistability is intimately linked with the zero-crossing of the pump detuning, as it changes due to the blue-shift 
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contributions. For negative detunings, the shift pulls the polariton feature towards resonance with the pump, causing 
the population to grow super-linearly. However, when the detuning becomes positive, the polariton is pushed away 
from resonance. The pump polariton population must lie on this curve if only the pump mode is significantly occupied. 
However, if other modes have finite occupation, as in the OPO state, this restriction no longer applies. Hence the 
conditions under which the OPO switches on can be found by determining which parts of the pump-only curve are 
unstable. 

The stability of the pump-only state is determined by considering the spectrum of small excitations with frequency 
u) and wave-vector q, of the form 



_ I — ifqx— ojt) I ■ i(qx— u;*t) 

p + se ' + le ' 



(2) 



which describes signal and idler modes with amplitudes s, i, and wave-vector/energies {kg = kjj — q, tUs = uop — lo) 
and [ki — kp + q, LUi ^ ujp + uj). Expanding to first order in s and i, this gives 

{ujo{ks) — ijs — ojp + oj)s + 2|pps + p'^i* — (3a) 



\Xs\' 
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(ujoiki) — i7i — i^p ~ uj*)i + 2\p\ i + p" s* = (3b) 



Considering these equations as an eigenproblem for the amplitudes .s and i* , the complex eigenvalues, w, are given 
by the condition that the determinant of the coefficients is zero, that is 

(u;o(fcs)-i7s -t^P+^^ + 2|X,p|p|2) X (4) 
( ujaih) + i7, - - CO + 2|X,np|2) = |x,p|X,np|4 

The threshold for instability corresponds to Im{cj} = 0, which occurs when 

A^ ^ —\Xs\^\X,\''\p\\ (5) 

Isli 

where F = 7s + 7,, and 

A = cjo(fcs) + uj^ih) - 2ujp + 2(|X,p + |X,n|p|2 (6) 

is the mismatch from the resonance condition 2cijp = ujf^ikg) + LJo(fci), modified by the blue-shift. Eqs. Q and ® 
provide a quadratic condition for |pp at the boundary of the instability region. 

Figdc) shows the instability region, plotted as a function of the signal angle 9s- For 6s — 6p, the extrema of 
the instability region correspond to the two turning points of the S-curve, so this analysis correctly predicts that the 
negative-slope portion of the curve is unstable. However, for the pumping conditions shown, the total extent of the 
instability region for all 0s is significantly greater: converting the extremal |pp values to pump powers gives a range 
between Ii and I2 on Fig^b), for which the pump-only state is unstable. 

The discussion above provides a criterion for determining whether the OPO state turns on for a given pump 
angle and energy: if there is a range of pump powers Ip for which the pump-only state is unstable, something more 
complicated, presumably the OPO, must turn on in between. In terms of Fig^b), this means that I2 > h- If, on 
the other hand, I2 < Ii, the system can jump straight from the lower branch of the S-curve to a stable state on the 
upper branch, and a simple bistability occurs. Of course, this does not rule out the possibility of the OPO switching 
on, but numerical experiments similar to those described in Ref.|0) only produce bistabilityiiSi A further possibility 
is that there is no unstable region, and |pp simply rises smoothly with increasing pump-power. Note that the OPO 
state does not require a bistable S-curve, as the unstable region of Fig^c) does not have to include 9s = Op. 

In Fig[21 the sectors of {9p, ujp) corresponding to these three types of behaviour are delineated. The figure shows 
that there is a minimum pump angle, ^ 10°, below which it is not possible to turn the OPO on. In the OPO sector, 
the threshold rises monotonically with increasing 6p and ujp, and there is no special significance to the 'magic angle', 
9p ~ 16°, at which the pump, signal and idler can all be made resonant with the lower branch dispersion. These 
conclusions are broadly consistent with the experimental results of Butte et a&. 

The key property which determines whether the OPO switches on appears to be the presence of a minimum in the 
higher threshold for 9s = 9p, as in Fig^c). If this is the case, the upper threshold, I2, is shifted above the higher 
knee in the S-curve, providing a gap in the stability of the pump-only curve. In Fig|31 it is shown that this correlates 
with the presence of a maximum in the 9s dependence of the bare mismatch when 9s = 9p: for 9p < 10°, instead 
there is a minimum for this value of 0,. 
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FIG. 2: Divison of the (6p,ujp) plane into sectors of different behaviour. In the 'stable' region, as the pump power, Ip rises, 
the pump population, increases smoothly. In the bistable region, there is a discontinuous jump in but no other state 
becomes occupied. The contours in the OPO region show the variation in the threshold, increasing in steps of 10(7^/^;) from 
the bottom left to the top right. The bare polariton dispersion is shown as a dashed line. 
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FIG. 3: The bare mismatch, ujoikg) + uooiki), plotted as a function of signal angle, 6s, for different values of pump angle. Op. 
The filled circles indicate the points where which 6a = 6p. For 6p < 12°, the dashed curves, these points are minima and there 
is no OPO. For 6p > 12°, they are maxima, and the OPO switches on. In this regime, the minimum mismatch always occurs 
close to 6s =0°. 



II. THE OPO STATES 



The previous section described the OPO states at threshold, when the signal and idler populations are still zero. In 
this section, the treatment is extended to include contributions of finite signal and idler populations to the dispersion 
renormalisation. This model is solved to obtain the variation of the pump, signal and idler populations above threshold 
for a range of signal angles. 

Eqs. ©J H3a|l and (|3b|) are only accurate to first order in s and i. For finite signal and idler populations, they are 
modified to include the renormalisation effects of the signal and idler populations, and the depiction of the pump: 
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-l-{u;o{kp) - i7p - iOp)p + + 2|s|2 + 2\t\')p + 2p*si = \^fp (7a) 
^ iLUoiks)-i-fs-i^s)s + i2\p\^ + \s\^ + 2\i\^)s + p^i* ^0 (7b) 
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{iuoih) - i7, - c^,)« + (2br + 2|sr + N )* + P = 0. (7c) 



Note the different form for the renormaUsation of each state: the self-interaction term has only half the strength 
of each cross term. In the OPO state the renormalisation does not produce a imiform blue-shift of the polariton 
dispersion. 

With the condition that uig and uJi = 2ujp — cug must be real, Eas. (|7bl7c() require 



^ _ r, . I'L \ 1 \v \2/n\„\2 , I „|2 , ol„-|2 

and for this value of uj. 



^[cuoih) + \X,\^i2\p\' + \s\' + 2\z\')] (8) 



(9) 



Then |sp and |ip satisfy Eq[Slwith the mismatch modified to include the signal and idler contributions to the 
renormalisation, 

A = c^o(fcs)+wo(fc0-2wp + |X,|2(2b|2 + |,5p + 2|in 

+ |X,p(2bp + 2|sp + |zn, (10) 

which, using Eq|51 represents a biquadratic equation connecting |pp and |sp. A further relationship between the two 
variables can be obtained by solving EqEc|for i in terms of s*, and substituting in the pump equation (|7a|l to get 

j^{ujoikp) - ijp - ujp)p + {\p\^ + 2|,s|2 + 2|z|2)p 

These two relationships can readily be solved to obtain the values of |pp and |sp in the OPO states, including the 
contributions of the signal and idler populations to the blueshifts. 

Such solutions for the OPO state are shown in Fig^ Above threshold, there is generally a region where the pump 
population, remains nearly flat, while the signal, |sp, increases with pump power. In a treatment with just the 
pump contribution to the renormalisation, |pp would be completeley flat for all powers above threshold. However, 
in the present model, a power is reached where the flat behaviour ceases, and |pp increases rapidly towards the 
pump-only line, at which point the OPO switches off, with |sp falling to zero. Examining the curves for different 
pump angles 9s reveals that the maximum Ip at the turn-over occurs when Og ~ 1°, and it decreases rapidly on either 
side, so at —7° and -1-9° there is only a small range of Ip in which the OPO state exists. 

The origin of this behaviour can be understood by considering the mismatch A, from Ea llOl which is zero when 
the pump and the renormalised signal and idler satisfy the triple resonance condition. As is shown in Fig^Jb), A is 
negative in the flat region, where the increasing blue-shift is thus pulling the OPO towards resonance, and helping 
the signal to increase. The turn-over occurs around the pump power where A becomes positive, so increases in the 
blue-shift drive the system away from resonance, and the OPO rapidly turns off. The dependence on signal angle 
is therefore a consequence of the different values of the bare mismatch 0^0(^3 ) + ^oiki) — 2ujp, which, because of the 
dispersion shape, always has its minimum close to 0^ = 0. Indeed, it can be seen from Fig|21 that this occurs at 
approximately -fl° for Op = 16°, and changes very little over the range of pump angles for which the OPO switches 
on. Thus, though the signal angle with the smallest bare-mismatch has a lower threshold, when the pump power is 
increased the growing blue-shift quickly switches that state off and the signal moves towards 9g — 0. 



III. OPO STABILITY 



The solutions for the OPO states described in the previous section are only of physical relevance if they are stable. 
This section describes a linear stability analysis of the OPO, and shows that there are stable solutions, but they 
become unstable when the pump population curves on Fig0|turn over. 



\ ^ 

50 100 

Pump Power, Ip ('] 



FIG. 4: The OPO states including signal and idler renormalisation contributions, (a) Relationship between the pump popu- 
lation, IpI^, and pump power, Ip, for signal angle = 0, and the corresponding signal population \sf . (b) The mismatch A 
from Ea llUI as a function of Ip for Os = 0. (c) The pump population variation, as in (a), for a number of values of Og. The 
curves for negative values of 9a are shown dashed for clarity. 



The stability analysis for the OPO proceeds in a similar manner to that for the pump-only state. The spectrum of 
small excitations about the OPO state is calculated using a linearised expansion of fluctuations in the pump, signal 
and idler modes. For an excitation wave- vector q and frequency w, these fluctuations satisfy a set of coupled equations, 
from which the dispersion is obtained using the condition that the determinant of the coefficients is zero. Thus 



dp p'^ + 2si 2{ps*+p*i) 

p*2_|_2s*i* dp 2p*s* 

2{p*s + pi*) 2ps ds 

2p*s* 2{ps*+p*i) s*2 

2{ps*+p*i) 2pi 2s*i 

2p*i* 2{p*s + pi*) p*^ + 2s*i* 



2ps 2{p*s + pi*) 2pi 
"8 + pi*) 2p*i* 2{ps*+p*i) 
2si* p^ + 2si 



ds 

+ 2si 
2si* 



P 



2s* i* 



2s*i* 



(12) 



where 
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{u;o{kp -q)~ i7p -Ljp + io) + 2{\pf + \sf + - \p\' 
(ujoikp +q)+ i7p -ujp-^) + 2(b|2 + |s|2 + |j|2) - 



(13) 



with equivalent expressions for dg, ds etc. The terms in the determinant are, in order of the columns, the coefficients 
multiplying the amplitudes of the fluctuations in p, p*, s, s*, i and i* . 

In Fig|Sl the eigenvalues calculated from Eq^] are plotted. Fig[5fa) shows the real parts of the roots for all the 
calculated eigenmodes. The calculated energies, w, are the excitations from the OPO state, which has three energies 
and wave- vectors {ojp,kp), [oji^ki) and (uji^ki). Thus, in a similar manner to the discussion concerning Eq|21 each 
excitation corresponds to six absolute energies and wave- vectors, (kp ±5, ujp ± lo), {ks ± (7, 0;^ ± u) and (fc^ ± q, 

±a;). These are the values at which the excited states can be created from, or transformed into an external photon, 
and it is these energies and wave-vectors which are plotted on the figure. The resulting dispersion curves, which are 
highly complicated, could be observed as parametric luminescence in the presence of the OPO state; indeed, analogous 
features have been observed in ultrafast pump-probe measurements^ 
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FIG. 5; Real part of dispersion curves for excitations from the OPO state with 9p = 16°, ujp = —27 and 6s =0°. (a) shows all 
the branches for the unstable state with Ip = 60.87'^/K, (b) the branches which are actually unstable, (c) shows the unstable 
branches for the nominally 'stable' state with Ip = 67.27'' /k, as discussed in the text. The solid line is the bare polariton 
dispersion and the large points are the {0,u}) locations of the pump, signal and idler. 

FigEJb) shows only the modes which are unstable, that is Im(ci;)> 0. For this pump power there arc two types 
of unstable excitation: small wave- vector modes, represented by the short, near horizontal, sections of dispersion 
passing through the pump, signal and idler, and much larger wave-vector modes which give rise to the longer arcs. 
For example, the long arc close to the idler is actually an excitation from the signal, with wave- vector g ~ fci — fcg, 
and it is paired with the arc at ^ —30°. 

FigEJc) corresponds to a slightly higher pump power than (b), where the small wave- vector instability has dis- 
appeared. However the large wave-vector modes remain, and indeed are always found in the present calculations, 
suggesting that the simple OPO state is never really stable. In fact, this is known to be true: in the experiments of 
Tartakovskii et a/ji^the OPO state is found to be accompanied by additional weak satellites, s' and i', at ks' = 2ks — kp 
and ki' = 2ki — kp, corresponding to the scattering processes p + p ^ s' + i', s + i — *■ s' + i' , p + s —* s' + i and 
p + i — > s + i' . Such satellites can also be seen in the numerical simulations of Ref.(|3). These particular states 
are selected because the corresponding q is equal to fcp — fcs, a resonance condition which is not accounted for in the 
derivation of Eq^| However, the comparison suggest that the instability associated with the large wave- vector modes 
is benign, and does not grow very large before being limited by non-linear terms. Distinguishing such instabilities 
from the more catastrophic variety is not possible within a linear stability analysis of the type performed here, which 
only shows that an infinitesimally small fluctuation in one of these modes will initially grow. 

Adopting the criterion that only the small wave- vector instabilities are important, Fig|Bl shows the portions of the 
OPO curves with negative 9s, in Fig^Jc), which are stable. It can be seen that for a particular pump power, the 
OPO state with the lowest pump population |pp is stable. This is consistent with the discussion in Ref.(|3), where 
it was argued that if a state with lower |pp were possible, the higher |pp state would be unstable. However, in the 
present treatment, the stable region for a particular angle extends beyond where it has lowest This is because 
stability is concerned with the response to infinitesimal fluctuations, and for a given pump power and signal angle, a 
small fluctuation in another mode may decay, where a large excursion into that mode would be stable. 
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FIG. 6: Stable regions of |ppversus 7p, for various signal angles 9s, according to the criteria discussed in the text. The solid 
lines represent the stable states. 

IV. DISCUSSION 

The results presented in this paper show that a proper treatment of the polariton renormalisation leads to a model 
which predicts the pump conditions under which the microcavity OPO can be observed, and explains the tendency 
for the signal emission to be close to the surface normal at high pump powers. 

The OPO threshold behaviour has been shown to be a consequence of how the mismatch in energy between the 
pump, signal and idler modes, A in EqlHI changes when their energies are renormalised by the self-interaction terms. 
This suggests an explanation for the observation, in ultrafast pump-probe measurements, that the 'magic-angle' for 
pumping is important: when an external probe is used the energies and angles of all three modes are determined by 
the experimental geometry, so the system cannot re-arrange itself in the same way. Hence, a strong response is only 
obtained when the bare polariton dispersion is resonant with the three modes, which gives the magic-angle condition. 
Indeed, as CW experiments are always carried out with the pump positively detuned from the bare dispersion, it is 
probably not surprising that the triple resonance condition is not important. 

The model also predicts that in the OPO region, there is a maximum pump power above which the OPO switches 
off, as well as a lower threshold. This behaviour is observed in numerical simulations such as Ref.|0), but not in 
experiments. This is probably a consequence of the inevitably inhomogeneous pump intensity associated with a finite 
excitation spot. A two dimensional numerical modeUfi with a Gaussian excitation profile shows that, as the power is 
turned up, the signal switches off at the centre of the spot, but there is always a region further out where the intensity 
is in the right range for the OPO to be active. 

A further prediction is that for given pumping conditions, there are a number of stable OPO states corresponding to 
different signal angles. This suggests that the actual state obtained will depend on exactly how the OPO is switched 
on, and possibly on random factors such as noise in the system. It may also be possible to steer the signal, by using 
a second probe beam to favour a particular emission angle. 
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